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Abstract 

The magnetic properties of matter under extreme conditions are of particular importance to un- 
derstanding the neutron star interior. One contributing factor to the magnetic field of a neutron 
star could be the ferromagnetic phase of nuclear matter. In this letter we present a self-consistent, 
relativistic description of ferromagnetism in dense matter, within which the ferromagnetic phase 
diagram for neutron matter is calculated. 
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The properties of dense nuclear matter are of interest to various branches of topical research. 
One area of interest is the description of the constituent matter of neutron stars. Neutron stars 
are extreme laboratories: extremes in pressure, density and magnetic field are just some of the 
factors that need to be considered when studying neutron star matter The neutron star interior is 
believed, at least in part, to consist of nuclear (baryonic) matter. Since charge neutrality of the 
star is required (a nett charge would result in a Coulomb force that will rip the star apart) the 
majority of the baryonic matter should be neutrons fl']. 

The origin of the very strong magnetic field of neutron stars has been a topic of discussion for 
some time. Brownell and Callaway ||2[] as well as Silverstein |3] proposed that the ferromagnetic 
phase of nuclear matter can make a significant contribution to the magnetic field of a neutron 
star This idea was investigated by various authors, using both non-relativistic and relativistic 
models, most recently Modarres and Pourmirjafari Ul, which also provides references to previ- 
ous investigations. The outcome of these investigations seems to be inconclusive regarding the 
existence of a ferromagnetic phase of neutron matter and appear to be very model-dependent. 

Due to the extreme densities encountered in the neutron star interior relativistic effects be- 
come more pronounced and thus it has been argued P, '5] that a relativistic description for the 
interior of the star may be necessary. Thus, to investigate ferromagnetic effects in the baryonic 
part of the neutron star interior a relativistic description of ferromagnetism in neutron and nu- 
clear matter needs to be considered. Such a description of ferromagnetism was investigated by 
Maruyama and Tatsumi [6], which based their work on that of Niembro et al (see |7] and ref- 
erences therein). We aim to further this work by presenting a self-consistent calculation of the 
magnetic field. We believe that a self-consistent approach is well-suited to this problem since the 
origin of the magnetic field in the ferromagnetic phase is the nucleon magnetic dipole moments, 
which in turn reacts to the presence of a magnetic field. 
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We perform this self-consistent calculation in the relativistic mean-field approximation. This 
is done by coupling the magnetic dipole moment of the neutron to the magnetic field of the 
ferromagnetic phase within the relativistic description pioneered by Walecka [IS-] called Quantum 
Hadrodynamics (QHD). QHD is an effective field theory model for nuclear matter with mesons 
as degrees of freedom. It had undergone various extensions and modifications, as reviewed in 
iQl and has been used extensively to study the properties of nuclei, nuclear matter and neutron 



stars; see for instance Uj, |9|, [lOl [IJl . We employ the coupling between the dipole moment and 
the magnetic field which was introduced by Broderick et al [12.1 to investigate the interaction 
between the magnetic dipole moment of nucleons and an external magnetic field. Within the 
relativistic mean-field approximation we use the Euler-Lagrange equation to derive an equation 
of motion for the ferromagnetic magnetic field (|4|, which we then solve self-consistently. As far 
as we can establish such a self-consistent calculation using this coupling has not been performed 
previously. 

To describe the interaction between the neutrons we couple the nucleon field to a scalar 
(sigma) and vector (omega) of meson field and to the photon field tensor F'^^. The meson 
describes the long-range attraction of the NN-potential, while the lo mesons accounts for the 
short-range repulsion liSJ. This description is valid for cold, neutral neutron matter, since the 
neutron matter is considered to be at zero temperature and the included mesons have no charge. 
The Lagrangian is 

-C = -Coirac + -CkG + -Cproca + -CeM - gv^JnOft// - gh^F''''o-ftyifr + (1) 

Here Xoirac is the free field component of the Lagrangian for the neutrons, Xkg is the Klein- 
Gordon Lagrangian for the scalar mesons and likewise -Cproca is the Proca Lagrangian for the 
massive vector mesons. -Cem is the free-field component of the electromagnetic Lagrangian. 
The coupling of the mesons to the different nucleon densities is the standard coupling as in 
cr^y are the generators of the Lorentz group. If A" = there is no electric field and F'^^cr^y 
reduces to -2B,2, in a co-moving frame of reference, with B, the components of the magnetic 
field and E, the three spatial four-component spin matrices. 

For the coupling strengths gy and gs the values of the QHDl parameter set given in jsl] 
were used. These coupling constants are chosen to fit various properties of nuclear matter at the 
nuclear saturation point ( The strength of the nucleon and electromagnetic coupling 

gh is left as a free parameter with the units of the nuclear magnetic dipole moment iC fm). 

The magnetic field is solved for in a co-moving frame of reference chosen, without loss of 
generality, so that the magnetic field lies in the z-direction. 

In the relativistic mean-field (RMF) approximation the nucleon and meson operators are 
replaced by their groundstate expectation values |5]. The equations of motion of the various 
fields are then found to be 
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For the equation of motion of the co meson field in Q, the RMF approximation results in the 
spatial components being zero and thus only the zeroth component is considered lUt] . 
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From Q it can be seen that the equation of motion of the nucleons is a modified version of 
the Dirac equation: 



[id, - g,m\lJ = [ - /or ■ V + gbBJil., +I3m*\lj, 



(6) 



with m* - m - gs(f> the effective mass and or and p the Dirac matrices. 

Thus the nucleon (Dirac) field (/r can be constructed in analogy to the free particle solutions 
of the Dirac equation ISJ. Therefore t/r is assumed to have the form 



ik-x-ie(k,/\)t 



(7) 



with e{k. A) the single particle energy. Since a non-zero magnetic field would break the spherical 
symmetry of the ground state it is convenient to express k in cylindrical coordinates as 



k^{k^,k,)^[^kl+kl, k,). 



(8) 



As pointed out by previous authors, l|6t] and references therein, care has to be taken when consid- 
ering the magnetic dipole moment (spin) of nucleons in the context of the Dirac equation. This 
is due to the fact that the spin operator does not commute with the Dirac Hamiltonian, the matrix 
on the right hand side of and therefore spin is not a good quantum number Using the form 
of the wavefunction in (|7]i, the single particle energies are found to be 



e(k, A) - g,MQ 




(9) 



Here e{k,A) is labeled hy A - +1, which refers to the contribution of the magnetic field to the 
single particle energies. A will serve to differentiate between the two species of neutrons, but this 
should not be mistaken for a spin label. As in the case of the unmodified Dirac equation both 
positive (particle) and negative (antip article) energies are found. 

Using the notation of Q the unnormalised Dirac spinors for particle are 



iff(k, A) 



k]^+[e(k,A)+g,,B^+m']lm' +A ■\^k\+m'^] 
[k:,+iky][e{k,X)+g},B;+A ■sJk-^+m'-] 



e(k,X)+g,,B-+A ■\lk\+m 
k~^[m'+A ^k\+m'-] 



{k^+iky\{e(k,X)+gbB^+X -^k^+m'^] 



1 



(10) 



The spinors are normalised so that 

A ip[, ^,{x) i/^k,A(x) = 6(k' - k)6A,A'- 



(11) 



By inspection of the single particle energies (|9]l it clear that when the magnetic field is zero, 
the normal dispersion relation is recovered and the lowest single particle energy state that a 
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Figure 1 : Illustration of the single particle energy levels that neutrons can populated in the presence of a magnetic field. 
The unpaired neutrons are of course the source of the magnetic field. 

fermion can occupy is the effective mass, m*. In the presence of a magnetic field the lowest 
single particle energy would be m* + ghB, for A - +1. The magnetic field therefore lifts the 
degeneracy between the states of the two neutron species. This splitting is illustrated in Fig. ^ 
for positive energy states, but also applies to the negative energy states. However the groundstate 
we consider only has positive energy nucleon states, which is filled to a certain energy, the Fermi 
energy. 

Using (|7]) to calculate the nucleon densities on the right hand side of equations (|2]-|4|, it can 
be shown that these equations are equivalent to the minimization of the energy density of the 
system, e. 



with jt/ the chemical potential and a step function to ensure that only energies below the Fermi 
energy are considered. 

We used the equations of motion to calculate the values of the different fields at a given 
density and value of the coupling gi,. The coupling constant was increased at a fixed density 
until a non-zero value of the magnetic field was found. The results are shown in Fig. In this 
figure the solid line represents the coupling strength gt, at a specific density which would effect a 
ferromagnetic phase transition when no mesons are present. The dashed line represents the case 
when mesons are included. 

We note that the vector mesons do not influence the phase transition since they couple to the 
nucleon density and merely shift the energy spectrum. This can be seen in the expression for the 
single particle energies 

The effect of the scalar meson is to reduce the mass of nucleon resulting in an effective mass 
m* . This reduction in the mass implies that for the same density the coupling constant gb has to 
be stronger than in the case when no mesons are present. Fig. (|2| clearly shows that the inclusion 
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Figure 2: Ferromagnetic phase diagram in density and coupling constant, gi,. The solid line represents the phase boundary 
when no mesons are included in the description and thus the neutrons only interact with the magnetic field. The dotted 
line represents the phase boundary when mesons are also included. The X indicates the point where gi, is at the value of 
the magnetic dipole moment of the neutron at a density which is equal to nuclear saturation density po = 0. 16 fm^^ . 

of the scalar mesons result in the phase transition occurring consistently at larger values of gi,- 

The shape of the curve in Fig. (|2]l at low densities and with mesons included, results from the 
interplay between the contributions from the scalar and nucleon fields to the energy density. 

From Fig. (|2]) it is apparent that for the ferromagnetic phase to present itself in cold neutron 
matter the coupling between the nucleon field and the magnetic field has to have some den- 
sity dependence. The coupling at normal nuclear density (the magnetic dipole moment of the 
neutron) is just not strong enough to affect a phase transition even at very high densities. The 
density dependence of gh might be calculated by including renormaUzation effects, but this was 
not investigated in the current work. 

It is possible that the densities at which the ferromagnetic phase of neutron matter is stable 
can be reached in neutron stars. In the calculation done here such a stable phase generates a mag- 
netic field of the order of 10'" Gauss. This magnitude of the field in itself is not enough to explain 
the origin of the magnetic field of neutron stars, which are calculated to have surface magnetic 
fields of the order of 10" - 10'^ G fis"]. The model presented here is a first approximation for 
neutron star matter, since only neutrons are included and only infinite matter considered. Further 
work would entail studying the effects of including other hadrons and leptons in the model to ob- 
tain a more realistic description for neutron star matter With such a model it might be possible 
to investigate the effects of a ferromagnetic phase on the equation of state of neutron star matter. 
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